Temporally stable coherent states in energy degenerate systems: The hydrogen atom 
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Klauder's recent generalization of the harmonic oscillator coherent states [J. Phys. A 29, L293 
(1996)] is applicable only in non-degenerate systems, requiring some additional structure if applied 
to systems with degeneracies. The author suggests how this structure could be added, and applies 
the complete method to the hydrogen atom problem. To illustrate how a certain degree of freedom 
in the construction may be exercised, states are constructed which are initially localized and evolve 
semi-classically, and whose long time evolution exhibits "fractional revivals." 
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I. INTRODUCTION 



Since their early introduction to quantum mechanics, the harmonic oscillator coherent states jlj have served many 
purposes. In Schrodinger's conception [Q, they were viewed as quasi-classical objects, purely quantum in definition 
though remarkably classical in behaviour. From this perspective, some authors have used these states with a classical 
limit to study correspondence between quantum and classical perturbation series ||. Stemming from Glauber's study 
t— ( ' j| , coherent states have a wide application in quantum optics [|j , primarily in representations of the electromagnetic 
field. In fact, being continuously parameterized, coherent states figure prominently in the theory of continuous 
representations p,0|. 

Many generalizations of coherent states appear in the literature. Each generalization tends to preserve a small 
number of the properties of the harmonic oscillator coherent states in the general scheme at the expense of the 
remaining properties. A recent generalization due to Klauder Q preserves many, at the expense of few. Klauder's 
-H i generalization gives states which (a) evolve among themselves in time (temporally stable), (b) are continuously 
parameterized, and (c) admit a resolution of the identity. As such, no reservations are made for "semi-classical" 
properties such as minimum uncertainty, though a certain degree of freedom to be discussed below remains within 
the construction which may be optimized according to additional concerns. Two studies have since appeared [p|jic|l 



proposing fourth conditions which eliminate this degree of freedom to be discussed in Section IV. 

As initially presented, Klauder's construction is appropriate for systems without energy degeneracies. With no 
additional structure, the resolution of the identity fails for degenerate systems. Energy degeneracies arise when 
independent operators commute with the Hamiltonian, suggesting a Lie algebraic approach to impose the additional 
structure. Thus in the presence of degeneracies, excepting those few cases of truly "accidental" degeneracies, the 
, Perelomov approach to constructing coherent states for the degeneracy group is an obvious and general path to take 

I' 

The attempt to generalize harmonic oscillator coherent states for the hydrogen atom problem is not new. Many 
involve the construction of the complete Perelomov states via the dynamical group SO(4,2) p^-[Tq|. Others, including 
the current approach, make some use of SO (4) coherent states [ ^6| , pT| . Others involve the construction of "temporally 
stable" coherent states: Klauder's original paper ||, a study by Majumdar and Sharatchandra |^], and another by 
Fox [ fl8| . The current construction differs from these in some key respects which will be pointed out in the course of 
the paper. 

Section || discusses the construction due to Klauder, and develops the extension of this construction to degenerate 
states. In Section [II, the construction is applied to the hydrogen atom problem. Section [fyjis a discussion regarding 



role that these states and other generalized coherent states hold in physical theory, namely, as representations rather 
than as physical states. The dynamics of an individual state are explored in Section [y], in which hydrogenic states 
are constructed which exhibit fractional and full revivals, and the main points of the paper are summarized in Section 



VI 



II. THE CONSTRUCTION 



The generalization due to Klauder ||] is applicable to the discrete portion of the spectrum of a Hamiltonian H. 
(A continuum generalization has also appeared fiofl .) So, for a non-degenerate H with eigenstates \n) and eigenstate 
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energies e n , the coherent states are given by (using atomic units) 

\ an) = M ^)± ^tr^ \n), (2.1) 

n VP™ 
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where s > and 7 is real. The factors p n are the moments of a function p(u) > 0, u > 0, 

/» CO 

p„ = / u n p(u)du, (2.2) 
Jo 

and the normalizing function M(s 2 ) is chosen such that 

00 2n 

l = ( s , 7 | s , 7 )=Af 2 (,s 2 )^— . (2.3) 

If the discrete portion of the spectrum is finite, the upper limits of the above sums (and the appropriate expressions 
hereafter) may be replaced with n max . 

The choice of p(u) is the remaining degr ee o f freedom up to the following restrictions: All the moments p n exist 
(up to n ma x if finite), and the sum in Eq. (2.3) exists for all s > 0. Such functions p(u) are known to exist: In the 
harmonic oscillator, p(u) — e~ u leads to the standard harmonic oscillator coherent states. 

The states \s, 7) clearly are continuously parameterized. The states exhibit temporal stability: exp(— iHt)\s, 7) = 
|s,7 + t). Further, given p{u) and M 2 (u), let k(u) be defined by 

k{u)M 2 (u) = p{u). (2.4) 

Then, the coherent states satisfy the resolution of the identity, 

1 = f d/x(s,7)| S ,7){s,7|, (2.5) 

in which the integration is given by 

J dp{s, 7 ) = Hm ^ ds 2 k(s 2 ) £ dry. (2.6) 



The limit is necessary to accommodate possible incommensurabilities of energy levels. In proving Eq. ( |2.5| ), integration 
over 7 yields the Kronecker delta £> e „e mJ which may be identified with S nm only in the non-degenerate case. Also, 



the identity of Eq. (2.5) should be regarded as a projection operator onto the states contributing to |s, 7 ), i.e. the 
discrete portion of the spectrum. These properties in concert make these states most useful in the representation of 
arbitrary, bound, time evolved states. 

To extend the construction to energy degenerate states, replace Eq. fl2.l| ) with 

\s,7,x) 

s n exp(-z 7 e„) 

71=0 



N{s)2_^ — V d «ka;), (2.7) 
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where d n is the degeneracy of the nth energy level, \n, x) are the Perelomov coherent states |^l| for the degeneracy 
group G, and the normalizing factor N(s 2 ) is given by 

In each energy degenerate subspace of the Hilbert space, the Perelomov coherent states satisfy the resolution of the 
identity 

i w = d n vol(H) I dr)(x)\n,x)(n,x\, (2.9) 
J x 
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in which H is the isotropy subgroup relative to the fiducial vector in the construction of the Perelomov coherent 
states, X = G/H is the quotient space formed by the degeneracy group with the isot ropy subgroup, and the measure 



drj is induced from the Haar measure on the degeneracy group. The states Eq. ( |2.7| ) therefore satisfy the resolution 
of the identity 

i = y dfi(s,j,x)\s,i,x){s,i,x\, (2.10) 

with 



dn{s, 1 ,x)= (2.11) 
lim — / ds 2 k(s 2 ) [ drfvol(H) [ drj(x). 
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Since the states \n,x) are formed by superpositions over states which share a common energy eigenvalue e n , they 
are also eigenstates of the Hamiltonian and so evolve simply in time. Accordingly, the states |s,7, x) preserve the 
temporal stability property of the non-degenerate construction. 

Majumdar and Sharatchandra's construction |J also makes explicit use of the Perelomov construction of coherent 
states for the degeneracy group, and Klauder's construction |8) is less explicit in this regard. A significant difference 
between this and these other extensions of Klauder's construction is the treatment of the factor d n . Among other 
problems p9j Majumdar and Sharatchandra incorporate d n into the measure after the summation of the state, an 
operation of questionable justifiability. Klauder, using an adaptation of SO (3) coherent states, incorporates d n into 
the states after, and therefore affecting, normalization. Fox [181 constructs temporally stable coherent states, but 
with a Gaussian taking the place of s n exp(— i"fe n ) in Eq. (J2 - 1 ) - As such, a similar adaptation of his states may be 
effected by including a factor of \fd^ as in Eq. (p?Z|). 



III. SPECIAL CASE: THE HYDROGEN ATOM 



The group theoretical treatment of the hydrogen atom is standard in the literature p0|-p2[|. For the hydrogen atom 
problem, there are two realizations of the degeneracy group SO (4). One uses the elements of the angular momentum 
vector, Lj, and a scaled quantum Runge-Lenz vector, Aj, as generators of the group, whereas the other decouples 
these six generators into two sets, Mj = \{Lj + Aj), and Nj — h(Lj — Aj). In the second representation, one finds 
that SO(4) = SO(3) ® SO(3), so that, loosely speaking, a Perelomov coherent state for SO(4) may be given by the 
direct product of two SO(3) coherent states. 



The SO(3) coherent states with the fiducial vector \j, —j), j — 0, h, 1, ■ 



are given by 
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The resolution of the identity for these states may be written 

2j + l f d 2 C 



1, 



(i + ICI 2 ) : 



-b\C)(j,CI, 



(3.2) 



where c£ 2 £ = dSRCrfSCi an d integration is over the entire complex plane projected from the unit sphere with the 
transformation ( = — tan | e~ 4< ^. In the hydrogenic realization, the representations of each copy of SO(3) are of equal 
dimension (n = 2j + 1) so the dimensions of the relevant representations of SO(4) are n 2 , n — 1, 2, 3, . . .. 

It is now straight-forward to construct the coherent states for the full system. The coherent states for the hydrogen 
atom problem by this construction are given by 



\s, 1X1X2) 

OO 

n(s 2 )y: 

71=0 



(3.3) 



s"exp(-z7e„ + i)(n + 1) 
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(2i)!Cf +mi d +m2 U,m 1 )|j,m 2 ) 



mi ,m,2 = —J 



l(j + m x )\(j - mi)!(j + m 2 )!(j - m 2 )!]V2(l + | Cl |2)i(l + ^pjj ' 



(3.4) 



The states | j, mi) |j, m 2 ) may be related to the standard hydrogen Hamiltonian eigenstates |n, I, m) via Clebsch-Gordon 
coefficients. The states Eq. (3.4) satisfy the resolution of the identity 



1b = -o / d ^{sn) 



d 2 Cid 2 (2 



(i + ICil 2 ) 2 (i + IC 2 | 2 ) 2 



|s,7,Ci,C2><s,7,Ci,C2|, 



(3.5) 



where the subscripted B is included to emphasize that this is more appropriately regarded as a projection operator 
into the bound portion of the Hilbert space. In the specific example of p{u) = e _ ", with moments p n = n!, explicit 
form may be given to N(s 2 ) and k(u) by 



2 \ - v-* 2 / 2 



(l + Zs' + s 4 )- 1 / 2 



and 



k(u) = 1 + 3u + u 2 



(3.6) 
(3.7) 



IV. SOME CLARIFICATION 



At this point, a few observations are in order. Primarily, the term "temporal stability" in no way refers to the 
time evolution of the structure in configuration space. Only through a rather generous interpretation does this 
construction "positively" solve the long standing problem of forming non-dispersing wave packets for the hydrogen 
atom. Temporal stability refers strictly to the mathematical property that the states evolve in time among themselves. 
With this property in mind, some authors Q have grossly overstated the the nature of the configuration space time 
evolution, while other authors |23 1^] have studied in detail the long-time evolution of individual states, even though 
there is no underlying physical basis either to provide for spatial coherence, or to presume states of this description 
are found in the laboratory at all. The question of how to prepare these states in the laboratory remains very much 
open. 

Much of the study of generalized coherent states rests more in the mathematical than the physical nature of 
mathematical physics. Glauber's motivation in his study of coherent states |Q| was not so much that coherent states 
are found in the laboratory, but that they provide a representation in which otherwise difficult calculations become 
feasible. Glauber noted that certain electric field operators have representations as sums over the modal annihilation 
operators. In diagonalizing these operators, one arrives at eigenstates of the modal annihilation operators. Restricted 
to a single mode, this corresponds to the annihilation operator definition of harmonic oscillator coherent states. Hence 
as annihilation operator coherent states, they arise from a representation, a point of mathematical convenience, not as 
a conclusion from the physics of the problem. In any case, generalizations of annihilation operator coherent states have 
appeared widely, though the physical motivation to study such definitions in any context besides as representations 
is unclear. 

Glauber also showed how these states may be constructed through the action of a displacement operator on the 
ground state. This definition was generalized by Perelomov [ p4| , a generalization which has been widely successful. 
This success is founded upon the properties of the dynamical group coming through into the set of coherent states, 
not from an assertion (which few researchers make) that an individual state by such a construction matches a state 
by some preparation in the laboratory. This success is of a mathematical, not physical, nature again resting upon the 
use of these states as a representation. 

Of Glauber's original three definitions, the approach which appears to invest the most physics is the minimum 
uncertainty construction. Indeed, squeezed states, a generalization of this construction, are used as descriptions of 
physical aspects of certain quantum optical experiments in the laboratory. Nieto et al. pf| have also developed a 
generalization which minimizes the uncertainty product of a pair of "natural" operators. Ehrenfest's relations then lead 
to the initial evolution of the quantum expectation values approximating classical evolution. Though this approach 
is strongest in terms of an underlying physical motivation, these states still lack (in general) any physical hypothesis 
which selects for states of this description in the laboratory. As an aside, the Nieto construction, though seen from 
time to time, is not as widely used as the Perelomov construction for perhaps two reasons. Firstly, though it is 
generally applicable in principle, many systems are intractable to carry through to completion (when the Hamiltonian 
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enters into the "natural" operators). Secondly, a certain tradeoff appears to be at work: This approach is somewhat 
less mathematically endowed than Perelomov's approach. 

Now consider Klauder's construction. All of the attractions are mathematical in nature. As initially presented, 
no reservations are made for coherence in configuration space (i.e. semi-classical behaviour) and there is no general 
physical mechanism which would result in finding these states in the laboratory. However, a certain degree of freedom 
remains in the construction, and two suggestions have separately appeared that a fourth requirement will simultane- 
ously eliminate the degree of freedom and ensure for the behaviour in configuration space. It is likely that a fourth 
requirement, if it exists, will be physical in nature. The requirement postulated by Majumdar and Sharatchandra || 
is that the measure found in the resolution of the identity corresponds to the "canonical" measure on classical phase 
space. They further assert that the measure uniquely identifies the set of coherent states. This assertion is false, as 
shown by Sixdeniers et al. [ p6| who demonstrate multiple measures corresponding to the same set of coherent states. 
Also, though it is convenient from a mathematical point of view, it is unclear why the measures should correspond 
at all from a physical point of view, or even if a meaningful identification (one to one) can always be made between 
individual coherent states and points in classical phase space. 

A fourth requirement is also postulated by Gazeau and Klauder [jiof which is motivated by an attempt to formalize 
the connection between the quantum parameters to the coherent state and the classical action-angle variables. Unfor- 
tunately, their requirement results in an angle variable whose rate of change with time is independent of the action, 
a rather special circumstance restricted to the harmonic oscillator and a small number of other systems. This is a 
severe limitation in terms of semi-classical behaviour, since this is clearly at odds with how the angle variable evolves 
in, say, the hydrogen atom problem. 

A degree of freedom also remains in Fox's construction [Q of Gaussian generalized coherent states, namely the 
width of the Gaussian in question. Note that in this context, the distribution in energy level, not configuration space, 
is Gaussian. Fox does not give any criteria which are intended to specify a suitable width. As with the Klauder's 
construction, this degree of freedom may be optimized according to the aim in mind. 

Hence in the absence of an acceptable fourth criterion (none is herein proposed), we carry on. This limits the 
construction to a mathematical tool, though an interesting mathematical tool it is. Note that the time dependent 
Schrodinger equation for a time-independent Hamiltonian in this coherent state representation becomes (in atomic 
units) 

|-( s , 7 |^) = -|< s , 7 |^> (4.1) 

whose solution is trivial. With this expression representing the unperturbed solution, this would make for an interesting 
point of departure towards a study of time dependent perturbation theory. That is, coherent states are useful when 
considered as an ensemble of states (a representation^, not as individual states. 



Speaking now in the specific case, some authors 23 19[] have suggested that the temporally stable construction of 
coherent states does not support the possibility of exhibiting full or fractional revivals as described by Averbukh and 
Perelman |27]] or Nauenberg . Firstly, in order to decide whether a state is to be found in the laboratory, a physical 
mechanism for the preparation of these states must be postulated. Until this has been done, the presence or absence 
of a phenomenon which is, after all, universal is not relevant. Secondly, these authors did not exploit the degree of 
freedom which remains in the construction. Without supplying a physical motivation which would lead to finding 
these states in the laboratory, we shall see that by exploiting this degree of freedom wave functions may be formed 
by the present construction which exhibit the full panoply of revivals. 

V. DYNAMICS 

Having thus constructed the states emphasizing, among other things, time evolution, it is now interesting to consider 
the behaviour of the states as evolved in time. Other authors have defined hydrogen atom coherent states with a 
variety of constructions and with various reports of evolution in "fictitious" time , or evolution along circular [^8 29 



v; 

or Keplerian elliptical orbits with possible, eventual state revivals Jl6|jl7[ ]. Coherent states also may be constructed 
by the present recipe which travel along elliptical orbits and exhibit fractional revivals. 

According to Averbukh and Perelman |27| , fractional revivals are a universal phenomenon exhibited by wave 
functions provided third order corrections and higher do not contribute significantly to a polynomial approximation 
to the energy eigenvalues over contributing energy eigenstates. Expanding about n = n, the hydrogen atom energy 
levels are 

e n = ~^= (5-1) 
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1 1 , 3 . — \2 ^ / — \3 

— ^ttt + — — n) — — : -(n — n) + — (n — n) +••■. 

2n 2 n 3V ; 2n 4 ^ ; n 5 v ; 

The first two terms lead to phase angles equal to multiples of 1-k for times in the vicinity of 

t = T r = ^(n) 4 . (5.2) 

This signals a wave function revival at t = T r provided the cubic term is small at the edges of the distribution, i.e., 
47r(An) 3 <C 3(n). Note that T r = T rev /2 in Averbukh and Perelman's notation. 

At the heart of these coherent states lies the function p(u). With p(u) — e~ u , one obtains a distribution in 
eigenlevels characterized by 

2 s 4 + 5s 2 +4 

{U) = 3 S 4 + 3 S 2 + 1 (5 - 3) 



and 

(An) = s wTiWTT ' (5 - 4) 

so that, taking leading order behaviour, An ~ U (n). Substituting into the above necessary condition for a revival at 
t = T r gives in^y (n) /3 <C 1 which is only satisfied for states in the immediate vicinity of the trivial coherent state, 
the ground state, and certainly violated by states with high principle quantum numbers. 

However, the function p{u) is a "degree of freedom" in the construction, and may be chosen according to application 
in mind. Accordingly, consider instead p(u) = exp(— u a ) for some constant a > with a view to constructing wave 
packets which exhibit strong revivals. The moments of this function are 

Pn = -T(^±l). (5.5) 

a a 

Many of the expressions involved in subsequent calculations may be handled using properties of the functions of 
Mittag-LefHer |3C]], though they will be treated instead by comparisons to expressions following from p(u) = e~ u . In 
fact, using p(u) = exp(— u a ) results in a set of coherent states closely related to those described by Sixdeniers et al. 

& 

Expressions for (n) and (An) 2 may be approximated by recognizing the scalings necessary to map expressions with 
a = 1 onto those with general a: 

n + l -> (n+ l)/a, s^s a . (5.6) 
Hence, one obtains to leading order (large values of s will eventually be involved) 

(n)-as 2Q , An -as", (5.7) 



so that An — yj a(n) . Substituting this into the minimal condition for the first full revival gives Ana 3 / 2 ^ (n) /3 <C 1, 
which may be satisfied if a is chosen sufficiently small. Without discussing the effect of changing a any further, there 
will be a tradeoff between large and small a: Large a will introduce many significantly contributing energy levels for 
a given (n) yielding good spatial localization but weak or non-existent revivals, whereas small a yields strong revivals 
of poorly localized states. Note that in this construction a small width in n follows from an appropriate choice for 
p(u) whereas the same may be accomplished by Fox's construction p8| by simply specifying the width to be narrow. 

As a typical example, consider the state depicted in Figure 0. For this state, a = 1/32 and s — 2.209 x 10 59 . 
This results in a state centred at (n) = 160 with a width of An = y5 for which one expects a full revival at 
t = T r — 1.373 x 10 9 . The parameters Ci,2 to the SO(4) coherent state were chosen provide an elliptical orbit with 
eccentricity e = 0.385, the major axis parallel to the x-axis, and the angular momentum parallel to the z-axis. The 
vertical axes are amplitudes of the wave functions on the x-y plane, calculated at the times indicated on a square grid 
80000 units in width centred at the origin. 

The evolution of this state is as expected. Initially localized, the state evolves semi-classically. The wave function 
then spreads out but remains close to the ellipse. As the expected times for the various fractional revivals arrive, the 
state exhibits the expected revival including the full revival at t = T r , even though an examination of the minimal 
condition for the first revival gives the debatable result 0.29 <C 1. 
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Figure |^ depicts the autocorrelation function for the same state as above, exhibiting the typical pattern character- 
izing revivals (compare with Figure 2 of Parker and Stroud |3l|| ) . Figure [| demonstrates how a small a shrinks the 
width of the state in n. 

Commenting again on the assertion by some authors that "temporally stable" states for the hydrogen atom can 
not exhibit this phenomenon p3| , [T9| . Their calculations involved, in terms of the present paper, a = 1, hence wide 
distributions in n which exhibited no appreciable revivals over the time frames calculated. Wave function revivals are 
a universal phenomenon depending on the nature of the energy eigenlevel spacings, and in the case of the hydrogen 
atom, the width An. By exploiting the fact that one may choose p(u) to one's liking, An may be reduced such that 
the resultant states do revive. Further, one study 23 used values of s leading to (n) equal to 25 and 400. With 
(n) = 400, from Eq. (5.2), T r ~ 5.36 x 10 10 though their calculation only extends to t = 5 x 10 9 . For the other study 
p9| , values of (n) exhibited were 10 and 200. A full revival should be found using (n) = 200 at about t = 3.35 x 10 9 
though their calculation only extends to t — 5 x 10 s . 



VI. CONCLUSION 



This approach to defining temporally stable coherent states for the hydrogen atom problem is similar in spirit to 
those of Majumdar and Sharatchandra ||, Klauder ||], and Fox |l8| ], Majumdar and Sharatchandra make explicit 
use of the SO (4) degeneracy group, but treat the factor d n differently. Although Klauder's construction is specific 
to the hydrogen atom problem rather than the present general approach, a similar term appears in that construction 
which disrupts normalization. Fox uses a somewhat different construction which avoids the use of the moments of 
some function p(u). Even so, one may construct coherent states in degenerate systems related to Fox's construction 
as the present construction is rela ted to Klauder's original work. 

The general construction of Eq. ( |2.7|) provides states with many useful properties. They form a complete set of states 
(in the bound portion of the spectrum) and evolve in time among themselves. This makes them a clear candidate 
for use in representations of time evolved, bound states. Further, there is a freedom in their definition which stems 
from the choice of p(u). If p(u) is appropriately chosen, hydrogen atom coherent states may be defined which exhibit 
the full range of phenomena exhibited by other approaches: initial semi-classical behaviour, interference between the 
head and tail of the state as it disperses about the Keplerian ellipse, localization on the Keplerian ellipse and wave 
function revivals at predictable times. 

The salient difference between these coherent states and other constructions is the is the natural and explicit manner 
in which the energy degeneracies are treated herein, via Perelomov's group theoretical construct of generalized coherent 
states. 
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FIG. 1. Amplitude ol coherent state on the x-y plane for times (a) t = 0, (b) t = T r /5, (c) t — T r /4, (d) t — T r /3, (e) 
t = T r /2 and (f) t = T r . 
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FIG. 2. Autocorrelation function of the coherent state at time t. For this state, the revival time is T r — 1.373 x 10 9 . 
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FIG. 3. Distributions in energy level for a — 1 (x) and a — 1/32 (+). 
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